Dynamics of the Jaynes-Cummings and Rabi models: old wine in new bottles 
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By using a wave packet approach, this paper reviews the Jaynes-Cummings model with and 
without the rotating wave approximation in a non-standard way. This gives new insight, not only 
of the two models themself, but of the rotating wave approximation as well. Expressing the models 
by field quadrature operators, instead of the typically used boson ladder operators, wave packet 
simulations are presented. Several known phenomena of these systems, such as collapse-rivivals, 
Rabi oscillation, squeezing and entanglement, are reviewed and explained in this new picture, either 
in an adiabatic or diabatic frame. The harmonic shape of the potential curves that the wave packets 
evolve on and the existance of a level crossing make these results interesting in a broader sense than 
only for models in quantum optics, especially in atomic and molecular physics. 

PACS numbers: 



I. INTRODUCTION 

Two fundamental models of quantum mechanics, pre- 
sented in most introductory textbooks, are the two-level 
system and the harmonic oscillator. Combining these 
two into a bipartite system gives many interesting mod- 
els, where two of the more studied ones are the Jaynes- 
Cummings ( JC) model 0, 0] and the Rabi model Q . The 
JC model was introduced for describing the interaction 
between a two-level atom and a quantized electromag- 
netic (EM) field, while the Rabi model was introduced for 
NMR systems. In this paper we will use the terminology 
of the JC model, a two-level atom and a quantized EM 
field, which, of course, does not restrict the results to such 
systems. The JC Hamiltonian is obtained from the Rabi 
one by simply imposing the rotating wave approximation 
RWA 4]. There are, however, special cases where atomic 
selection rules make the JC model exact; the RWA terms 
naturally vanish [f| . In this approximation, exact analyt- 
ical solutions exist, and, in spite the simplicity of the JC 
model, the dynamics have turned out to be very rich and 
complex, describing several physical phenomena. Among 
these are; Rabi oscillations d, |(|, collapse-revivals Q, 
squeezing Q, atom-field entanglement @, non-classical 
states as Schrodinger cats [l(| and Fock states [ll| and 
anti-bunching [12j. The JC model originally thought of 
as a single atom, single field mode interaction, has shown 
to be applicable for several other types of systems, for ex- 
ample, trapped ions [HI], Cooper-pair boxes [3], "flux" 
qubits and Josephson-junctions pjj]. With the ex- 
perimental progress of some of the above mentioned sys- 
tems, the coupling between the systems may be made 
very large, and the RWA breaks down so that only the 
Rabi model describes the dynamics correctly. As the 
number of excitations, field plus atom, is a constant of 
motion in the JC model, there are only two physically rel- 
evant parameters; detuning A between the atomic tran- 
sition frequency f2 and the field mode frequency lu and 
atom-field coupling go. In the Rabi model, however, the 
number of excitations is not conserved and all three pa- 
rameters are of importance. In general, the RWA is justi- 



fied for small detunings A = £1 — uj and small ratios of the 
atom-field coupling divided by the atomic transition fre- 
quency, go/&- In atom-field cavity systems, this ratio is 
typically of the order 50/^ ~ 1CT 7 -10~ 6 , 17]. Recently, 
cavity systems with very strong couplings have been dis- 
cussed [l8| . The ratio may also become order of magni- 
tudes larger in solid state systems [l4[ , and the full Rabi 
Hamiltonian, including the virtual processes (also caller 
counter-rotating terms), must be considered. The neglec- 
tion of the counter-rotating terms may have interesting 
physical consequencies on nonlocality and causality [H^ ]. 
The effect of the RWA in other systems has been stud- 
ied, and to mention a few, trapping of atoms by light 
[20| . driven Rabi model l2l[. multi- level atoms and/or 
fiel ds [23 ] . the micromaser [25| semiclassical Rabi mod- 
els [26l ] . transition between vibrational states in a water 
molecu le [ 23| . dipolar molecules [24j, light-matter inter- 
action [27fl. electromagnetically and self induced trans- 
parency [28|, [29( and laser modified collisions [3fJ • 

Existence of analytical solutions for special cases of the 
Rabi model has been discussed [3l[. As no general sim- 
ple solutions are known, much work has been devoted to 
various analytical approximate or numerical approaches, 
such as; perturbation theory [H, [H, l34l 139 . 1361 \3l\ , 
path-inte gral s [38| . continued fractions l39ll. variational 
methods Eh coupled cluster methods [4J|, displaced 
oscillators [42L |43j], app roximated unitary transformed 
Rabi Hamiltonians [44j and numerical studies [H, [46|, 
1471 [H , |49| . The time-dependent Rabi Hamiltonian has 
been considered in [H(J, and a particular equivalence re- 
lation between the time-dependent JC and Rabi models 
was obtained. Interestingly, the Rabi model has shown 
to be chaotic [5l|, contrary to the JC model [52| ■ The 
Bargmann representation [53j , similar to our approach, 
has been considered for both the JC [54[ and the Rabi 
model [55} . 

Using a different approach than the standard ones we 
discuss and review some of the work that has been done 
previously on the two models. The method used in this 
paper is a numerical wave packet propagation, which has 
been briefly applied to the Rabi model before p5l . l46| . but 
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to the best of our knowledge not to the JC model. The 
system Hamiltonians are usually represented in the field 
creation and annihilation operators and a, but here we 
instead work with the field quadrature operators P and 
Q. As these obey the standard canonical commutation 
relations, the system is equivalent to a particle, with unit 
mass, moving in two coupled potentials, and the dynam- 
ics is given by the evolution of some initial wave packet. 
The idea of this paper is to give a deeper understanding 
of the JC and Rabi models by consider the wave packets 
evolving on these potential curves. Since the JC model is 
exactly solvable, it might seem unmotivated to use this 
numerical study. However, from the wave packet picture 
it is easy to understand some of the systems behaviours, 
and, additionally, due to the "graininess" of the cavity 
field most analytically obtained quantities predicted by 
the JC model are given as infinite sums lacking closed 
forms. As in related problems in molecular physics [5^ |. 
it is often convenient to investigate the system in the 
diabolic or adiabatic basis [HI, Hy, |47], H7[. The proper 
basis used, depends on the system parameters and ini- 
tial state. The critical evolution takes place close to level 
crossings between the two energy curves, and such a level 
crossing exists in both the JC and Rabi model. The two 
energy curves are given, in both models, by two displaced 
harmonic oscillators, and they are coupled by a constant 
coupling in both cases, while in the JC model there is 
an additional " momentum" P dependent coupling which 
gives rise to very different dynamics. The presence of the 
P dependent coupling makes the wave packet evolution 
of the JC model less intuitive and more complex. Com- 
paring the two models, combined with the knowledge of 
the validity regime of the RWA, one may interpret the 
consequences of this coupling term as well as other ef- 
fects. In the intermediate regime, when one is neither in 
the adiabatic nor the diabatic validity regimes, a more 
detailed analysis is needed to understand the full evolu- 
tion. This is only briefly mentioned in this paper. 

The paper proceeds as the following. In the next sec- 
tion [ED we introduce the JC and Rabi Hamiltonians in 
their regular forms and also give them in the conjugate 
variable picture. In subsection III CI the corresponding 
curve crossing problem is mentioned and how to approx- 
imate it as a Landau-Zener problem 58]. The following 
section IIIII presents a background of an adiabatic ap- 
proach, the two different bases, diabatic and adiabatic, 
are defined, and the approach is applied to both the JC 
and Rabi model. The validity of the adiabatic method is 
numerically studied in Sec. IIVI using the split operator 
procedure to calculate the fidelity. In Sec. [V]we present 
the numerical results of our studies. Several plots are 
shown to give a deeper insight of the dynamics and the 
phenomena are discussed from a wave packet point of 
view. Known interesting effects, such as entanglement, 
collapse-revivals, squeezing and Rabi oscillations, are es- 
pecially studied. Finally we conclude the paper in sec. 
IVII with a summery. 



II. THE JAYNES-CUMMINGS AND RABI 
MODELS 

A. Introducing the JC model with and without 
RWA 

The simplest fully quantum mechanical model de- 
scribing light-matter interaction, naturally considers one 
atom, with a single atomic transition, interacting with 
only one mode of the field. Such an approach leads to 
the JC or Rabi models. The single atom transition and 
single mode approximation relies on small or vanishing 
coupling elements between other levels [E{|. Thus, math- 
ematically one is left with a two- level atom (spin- 1/2 par- 
ticle) coupled to a harmonic oscillator. For a dipole in- 
teraction, Hi n t = d ■ E/h, where d and E are the atomic 
dipole moment and electric field respectively, a micro- 
scopic derivation [5§| gives the Rabi Hamiltonian 

HrcM = ftw (ala + +^Y<J z +fi9o (c + + o~) (a f + a) . 

(1) 

Here, (a) is the creation (annihilation) operator for 
the field mode; a\n) = y/n + l\n+ 1) (a\n) — y/n\n — 1)), 
the sigma operators are the standard Pauli matrices act- 
ing on the two-level atom; a z \±) = ±|±), o p± |±) = |=f)> 
a x = a + + <t~ and a y — —i (<j + — <j~), lu (fi) is the 
field (atomic transition) frequency and go is the effective 
atom-field coupling. In deriving |T]), the dipole approx- 
imation has been assumed, where the variation of the 
field is neglected on the atomic length scale, and the ki- 
netic energy of the atom is omitted, valid for large and 
moderate temperatures [60(. 

The interaction part contains four terms; a + a i (<r — a) 
simultaneous excitation (de-excitation) of the atom and 
field, and a + a (er~a^) excitation of the atom by ab- 
sorption of one photon (de-excitation of the atom by 
emission of one photon). In the Hciscnberg picture, 
the operators evolve (free field and free atom case) as: 
a(t) — aexp(icut), tz^(t) = exp(— iu>t) and a (t) = 
a ± exp(±ir2t). Thus, the interaction terms precess with 
either the frequencies |f2— u\ (energy conserving terms) or 
f2 + u> (non energy conserving/counter rotating terms). 
Most often |f2 — oj\ <C fi + u and the fast oscillating 
terms are rejected from the Hamiltonian, resulting in the 
Jaynes-Cummings Hamiltonian 

Hjc = hio fc^a + + ^-o- z + hg (a + a + cr~a t ) . (2) 

However, for a large atom-field detuing (but still within 
the atomic two- level and single mode approximations), 
the above constrain may be violated and the non energy 
conserving terms can not be excluded. As mentioned in 
the introduction, also the ratio go /O is important for the 
validity of the above RWA. The two conditions for the 
RWA must in general be met simultaneously, meaning 
that near resonant interaction does not gurantee validity 
of the RWA if not go < and vice versa. 
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B. The models in the conjugate variable picture 

The JC model is usually given and solved within the 
"number" -basis of the field {|n)|±)}. Here, however, we 
will work in the conjugate variable basis of the field. We 
introduce the conjugate variables P and Q related to the 
creation/annihilation operators as 



(3) 



In the transformed basis the Rabi ([T]) and the JC @ 
Hamiltonians read 



^(at-o), 



Q = \h W 
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Relating P and Q as momentum and position opera- 
tors respectively, the above Hamiltonians can be inter- 
preted as a two-level particle, with quantized position 
and momentum, confined in a harmonic trap and in- 
teracting with a "classical" field with mode "variations" 
g(Q,P) ~ Q or g(Q,P) ~ (Q + The Rabi Hamil- 
tonian could be obtained, for example, by considering a 
harmonically trapped two-level atom driven by a classical 
field, whose wave length is much longer than the extent 
of the atomic wave packet in the trap and with a node at 
the centre Q = 0. In other words, there is a direct con- 
nection between the above example of a trapped atom 
and an atom interacting with a quantized field. Relat- 
ing the JC Hamiltonian to similar systems is less trivial, 
since the atom-field coupling is momentum dependent. 



C. The corresponding curve crossing problem 

The structure of the Rabi Hamiltonian (0| is more 
easily seen by applying the unitary transformation U = 
{ a x + c z ), giving the transformed Hamiltonian 



HRabi — U'HjiabiU 



p2 



-Q 2 



goVZruvQ 

Ml ,-— 

— -goVZruvQ 



(6) 



By completing the squares we may write the Hamiltonian 

as 



ilRabi — -pr- 




%0 



(7) 



'/2 is the regular harmonic oscillator 



where Vh{x) 

potential. Thus, it is seen that the problem is equivalent 
to the one of a particle moving in two coupled equally, 
but opposite, displaced harmonic potential curves. The 
two energy curves cross for Q = 0, and the off diagonal 
coupling Ml/2 results in an avoided crossing. As is well 
known; the main population transfer between two levels 
occurs close to the crossing. Here the distance between 
the curves relative to the coupling amplitude is the small- 
est. Therefor, in some situations the coupled dynamics 
may be considered only in the vicinity of the crossing. 
By linearizing the model around Q — we get 



tl Rabi — 




(8) 



If the wave packet is well localized (the characteristic 
width of the wave packet Aq is small compared to vari- 
ations of the harmonic potential) , the " momentum" and 
"position" may be replaced by its classical counterparts 
at the crossing; P — > v and Q — > vT, where T is the time. 
Neglecting the constant terms, the Hamiltonian reduces 
to the one of the Landau-Zener model I5S 



PiRab 




(9) 



where the classical velocity v will, of course, depend on 
the initial state of the system. The time dependent prob- 
lem given by the Hamiltonian ([9]) is analytically solvable, 
with asymptotic solution for the population transfer to 
the other level given by 



Plz = 1 - exp 



ntl 2 



(10) 



The validity of this semiclassical model may be studied 
in a detailed way similar to the one in [6l| . 

The same procedure may be applied to the JC Hamil- 
tonian (O; 



Hjc — - 



P 2 




_M 
4w ' 

(11) 

One notes that the displaced oscillators are shifted half 
the amount as in the Rabi Hamiltonian, and that there 
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is an additional P-dependent coupling of the oscillators. 
As the RWA implies that non energy conserving terms 
are omitted, it follows that the number of excitations is 
conserved. Thus, N = ^ + ^Q 2 + \<j z commutes with 
Hjc, and one may work in a rotating frame with respect 
to N. Applying the same unitary operator to the new 
interaction picture Hamiltonian H' JC = Hjc — htoN, one 
finds 



H' 



JC 



,9o 



hw hA / h 



hA 



«5o 



p 



-go 



Q 



(12) 



The energy transfer between the atom and the field in 
the JC model is give by H' JC , and from (fi"2"| we note the 
similarities with the Landau-Zener problem but with the 
addition "momentum" dependent coupling. An interest- 
ing observation is that if the wave packet approaches the 
curve crossing with a high momentum, an adiabatic evo- 
lution is still possible due to the large curve couplings. 
At the one hand, "velocity" determines the steepness of 
the crossing curves, but on the other is also determines 
level separation close to the coupling. This peculiar curve 
crossing model is, of course, an artifact of the RWA, but 
none the less it may give insight into other related sys- 
tems, for example atomic and molecular scattering pro- 
cesses 1621. 



III. THE ADIABATIC APPROACH 
A. Review of the adiabatic principle 

The concept of adiabaticity is widely used in a va- 
riety of physical and chemical systems, while, what is 
known as the adiabatic theorem (63j is only stated for 
time-dependent Hamiltonian systems. Non the less, the 
principle may be transferred to other systems where the 
dynamics is governed by the evolution of some initial 
state (64[. Here we review the results of (65|, and apply 
it to our models in the following subsections. 

Given a Hamiltonian of the form 



H = 



P 

2m 



V(x) 



fiA(x,p) 
2 

hg(x,p) 



hg(x 1 p) 
HA(x )P ) 



(13) 



The diagonal elements are often refereed to as detuning, 
and for simplicity they will be taken constant A(x,p) = 
A, and the off diagonal ones are level couplings. We 
introduce the parameter 9 accordingly 



tan 29 = 



2g(x,p) 



such that the unitary transformation 



U = 



cos ( 
sin f 



— sin 9 
cos 9 



(14) 



(15) 



diagonalizes the last term of the Hamiltonian (fi"3| with 
corresponding eigenvalues 



(16) 




9 2 (x,p), 



and eigcnstates, defining the new adiabatic basis states, 



IT> = 



cos ( 
sin t 



— suit 
cos 9 



(17) 



However, since in general [p, g(x,p)] and 

[V(x), g(x,p)] ^ 0, the new Hamiltonian will not be di- 
agonal. Introducing the two functions / = cos 9 and 
h = sin#, and using the identity 



U^pU = p- a y hd9 : 



where 



89 = — 

ox 



the transformed Hamiltonian takes the form 

H = H a d + H cor , 
with the adiabatic part 



H ad = ^- + (hd9) 2 + fV(x)f+hV(x)h+ 
and the correction part 



A 

-A 



(18) 
(19) 
(20) 

(21) 



TT _ 1 

-"cor — o„, 



h 2 d 2 9 - 2ih{d9)p 

-H 2 d 2 9 + 2ih(d9)p 



+ {fV(x)h-hV(x)f)ia y . 

(22) 

Note that if V(x) commutes with g(x,p) we have 
fV(x)f + hV(x)h = V(x) and fV(x)h - hV(x)f = 0. 
As it will turn out, this will be the case in our models, 
and we therefor assume commutability. 

An adiabatic evolution is obtained if the dynamics is 
dominated by the diagonal adiabatic Hamiltonian H a d 
of (|2TJ)) . while H cor only slightly affects the propagation. 
Clearly, this depends on the smoothness of 9 and the 
momentum p, but also on the wave packet shape. Using 
the definition of the angle parameter (jT4j) we find 



89 = 



d 2 



Adg 



A(d 2 g (A 2 +4g 2 )-8g(dgf 



(23) 



(A 2 



V) 2 



Under the adiabatic approximation, the state is evolved 
on the adiabatic energy curves s± — V ± A according to 



* ad (x,i) 



exp 



exp 



1 \ 2m 



<pf(x,o)\ t) 



£_ t 



^fixM I)- 



(24) 
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In the opposite limit, the non-adiabatic case, the evolu- 
tion takes place on the diabolic energy curves (also called 
bare energies), defined by the diagonal elements of the 
Hamiltonian. Note, though, that the diabatic energies 
depend on the basis used. 



B. Application to the JC model 

Before proceeding, we introduce scaled dimensionless 
variables. We scale energies with the photon energy hu>, 
"lengths" Q by \Jh/ui and time by 



t = wT 



go 

.9o = — • 

U) 



(25) 



on the last two 



For convenience we drop the tildes 
scaled variables. 

Introduced in the subsection III C( it is convenient to 
work in an interaction picture with respect to excitation 
operator, 

A 1 
~2~ 



H' 



JC 



H 



JC 



N 



50 71 ( q ~ ip ) 



A 

~ ~2 



ip) 



(26) 

with the detuning A = — 1. As the interaction pic- 
ture Hamiltonian is lacking the kinetic energy operator 
and the " external potential" , the correction H cor to the 
adiabatic Hamiltonian vanishes. The eigenstates are 

| |)n = cos0Vn-l(?)|+) + sin 6tjj n (q)\-), 

(27) 

| 1)„ = sin0Vn-i(?)|+) - cos0ip n (q)\-), 

where 

2g oy / n 



tan 29 = 



(28) 



and tp n (q) is the nth eigenfunctions to the harmonic os- 
cillator. As q ± ip act as raising and lowering operators, 
the diagonaliced Hamiltonian in the original picture and 
within the nth block reads 




^ = £ + £ + ^.±1/1 = 
Thus, the eigenvalues are 

E±(n) = 



9o n - 




(29) 



(30) 



The "adiabatic limit" in the JC model is generally 
taken as A > goy/n- It has been shown [65[ that this 
limit governs adiabatic evolution, but that this is not the 
only possibility. Expanding the solutions (|2"T)) and (13"01) 
to first order in g y/n/A reproduce the familiar results 
from "adiabatic elimination" in the JC model (66j. 



C. Application to the Rabi model 

Applying the above method to the Rabi Hamiltonian 
(|4]) give the adiabatic Hamiltonian 



H ad ^ F - +V±(q), 



where 



V±(q) 



2 



2&gl 



(fi2 + 2 g$q^ 




2g 2 q 2 



(31) 



(32) 



and the corresponding adiabatic states given in (|17|) . In 
the limit Q — > 0, the adiabatic correction vanishes and 
the problem is analytically solvable, since (f3"Tj) becomes 
identical to two displaced disconnected harmonic oscilla- 
tors. 

The standard way of defining an adiabaticity criteria 
is by letting the "distance" between the adiabatic ener- 
gies become much larg er than the amplitude of the off- 
diagonal couplings [63J. In our particular model, we ex- 
plicitly have 



89 



d 2 9 



»ffoV2 
n 2 + 2g 2 q 2 ' 

4n g ^V2g 
(Q 2 + 2 5o V) 



(33) 



For smooth angles 9, the main non-adiabatic contribu- 
tion will come from the term 2(d9)p, where a large p is 
obtained by having a large (n) which is known to violate 
adiabaticity. We directly note that the standard adia- 
batic limit Q 3> go gives a small coupling. However, we 
also note that if q 2 ^> (v2<?o^ — fl 2 ) /2<7q, an adiabatic 
evolution is obtained, which is intuitive as we are far 
from the curve crossing. It is instructive to analyze the 
diabatic and adiabatic energy curves for different sets of 
parameters which is shown in fig. 1. Note that the dia- 
batic energies are just the displaced harmonic oscillators 
Vh(q ± V2go)- In (a), the atomic frequency ft is larger 
than the coupling and we have in this situation the reg- 
ular adiabatic limit. In (b)-(d) f2 is smaller or equal to 
go and an adiabatic evolution is only expected when the 
wave packet is well localized at positions \q\ ^s> 0. We 
note how the barrier in the middle at q = of the adi- 
abatic curves is lower, but wider, for increasing f2. This 
peak will most likely govern a non-adiabatic transition 
between the adiabatic states. This again, indicates that 
adiabatic evolution is only obtained far from the crossing 
in these situations. For large enough \q\ the adiabatic and 
diabatic states coincide for any parameters. When £l/go 
is small, a large \q\ is needed to assure adiabaticity, while 
a large ri/ffo induces adiabaticity for all "positions" q. 
Thus, since the states we propagate will most often pass 
through the level crossing we conclude that for > go 
the behaviour may be understood from the evolution on 
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the adiabatic energy curves, while for £1 < go it is better 
to use the diabtic curves. 

For simplicity, the initial state will be given by the 
atom in its excited state and the field in some photon 
distribution, hence 

0) = Ml, 0)|+) = MQ, 0) (cos0\ t) - sinfl| I)) . 

(34) 

This limits the analysis, since more effects may be seen 
with other initial conditions, but that would greatly ex- 
tend the leangth of the paper. For an initial Fock state, 
a^a\n) = n\n), or coherent state, a\v) = v\v), we have in 
g-representation respectively 



^n{q, 0) 



1 / 1 



\/2 n n\ V 71 " 



1/4 



H„(g)e , 



1 7T 



1/4 



(35) 



where H„ is the nth order Hermite polynomial. We note 
that for initial Fock states, adiabaticity is not expected 
if Q < go, since their wave functions are non-zero around 
g = 0. However, for odd n we have ip n {q = 0, 0) = 0, but 
numerics indicate that adiabaticity is lost also for such 
states. 
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Figure 1: Solid lines display the adiabatic energy curves (|32[) 
for four different sets of dimensionless parameters indicated in 
the plots. The dotted lines show the corresponding displaced 
harmonic potentials Vh(q ± v2ffo). Here the "momentum" 
p = l. 



the split operator procedure [67[, which relies on sepa- 
ration of the evolution operator into a momentum and 
position dependent one. The wave packet is either given 
in the q- or the p-reprcscntation, by a simple fast Fourier 
transform (FFT). The initial state of the field is given in 
eq. (IM)) . and it will be propagated by, either the adia- 
batic Hamiltonian (|3ip or the full Hamiltonian ([?]) . From 
the corresponding states, labelled ^ ad (q, t) and ^(q, i) re- 
spectively, we measure the validity of the approximation 
with the fidelity 



F(t) = 



$>*(q,t)$ ad (q,t) dq 



(36) 



Numerical simulations indicate that a large detuning, 
fl -C 1 or £1 ^> 1, governs an adiabatic evolution. These 
limits are opposite to the one related to validity of the 
RWA. However, the limit (go/iY) — > assures applica- 
bility of the adiabatic approximation and RWA. Non the 
less, there are parameter regimes where the adiabatic ap- 
proximation is justified, while RWA breaks. Particularly, 
for a fixed ratio go/Q, the limits Q — ► oo or Q — * guar- 
antee adiabaticity. 

In fig. [2] we show the dependence of the fidelity (1361) . 
for an initial empty cavity (Fock state), on the coupling 
amplitude go and time t. In the left plot il = 2 and in 
the right one $7 = 0.5, and it is clear that keeping go/^t 
small results in an adiabatic evolution. The next plot, 
fig [31 gives the same, but now for an initial coherent 
state of the cavity field with an amplitude v = 4. In 
(a) fi = 10 and in (b) £1 = 0.1. The colour scale is 
the same for both figures. By lowering the amplitude v 
the fidelity would be increased. It should be emphasised 
again that for a larger detuning between the mode and 
atomic transition frequency, one would also obtain a more 
adiabatic evolution. 
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IV. VALIDITY CHECK OF THE ADIABATIC 
APPROXIMATION 



Figure 2: (Colour online) Plots of the fidelity (|36p . for an ini- 
tial excited atom and empty cavity, as a function of scaled 
time t and atom-field coupling go. In (a) the atom transition 
frequency is Q = 2 and in (b) fl — 0.5. The colourbar to the 
right displays the value of F(t). 



The validity of the adiabatic approximation is studied 
from a full wave packet simulation. The method used is 




Figure 3: (Colour online) The same as fig. [5] but for an initial 
coherent state with (n) — \v\ 2 — 16. The atomic transition 
frequencies are fl = 10 and SI = 0.1. The colour scale is the 
same as in fig. [2] 



V. NUMERICAL INVESTIGATIONS 

In this section we investigate in a non-standard way 
some of the more well-known phenomena predicted by 
the JC model. Wave packet propagation will be pre- 
sented for both the JC and Rabi model. As mentioned 
in the introduction, one advantage of the method of using 
wave packet propagation is that it is valid for all param- 
eter regimes. However, one should keep in mind that 
the adiabatic approximation is only valid within certain 
parameter regimes as seen in the previous section. The 
numerical results will be analyzed in terms of the wave 
packets propagating on the appropriate adiabatic or di- 
abatic energy curves. 

A. Rabi oscillations 

Rabi oscillations Q describes the periodic population 
swapping between some few internal levels, driven by an 
external " field" . For the JC model, Rabi oscillations 
manifest themselves by population transfer between the 
upper and lower atomic state, equivalent to an energy 
transfer between the field and the atom. These oscilla- 
tions are " exact" when the field is initially in a Fock state 
and the systems are in resonance. Rabi oscillations have 
been measured in several JC type of systems [f| , and also 
theoretically studied for the Rabi model [68[. When the 
non energy conserving terms are included, the structure 
of the oscillations may be greatly affected. Figure H] gives 
two examples where circles correspond to the Rabi model 
and solid line to the JC one. In (a) we have ft < uJ,go, 
while in (b) > u>,go- The peaks of (a z ) in fig. [4] (a) 
come with a period of 2tt (in scaled units, where w = 1). 
The adiabatic frame is no longer appropriate for giving 
an intuition of the behaviour, but instead one may con- 
sider the wave packets as evolving on the two displaced 
oscillators. For the Rabi model, the wave packet starts at 
the origin and then splits up into two parts, respectively 
propagating down the two displaced oscillators. After 
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Figure 4: The figure shows the atomic inversion (Rabi oscil- 
lations) for an initial empty cavity and excited atom. Circles 
give the non-RWA results, while solid lines the RWA ones. In 
(a) fi = 0.2 and go — 2 and in (b) fi = 4 and go — 2. 



one classical period of oscillation the two packets return 
to their original positions and the interference gives rise 
to the peak. In the JC model, the wave packet will not 
split into two disconnected parts because of the p de- 
pendent coupling which tend to bound the wave packet 
to the origin. Therefor, there is no collapse in the Rabi 
oscillations. In (b), the situation is adiabatic and since 
the adiabatic energy curves, fig. [l](a), do not form two 
displaced oscillators, a collapse of the Rabi oscillations 
is not as clear as in (a) for the Rabi model. In fig. [5] 

(a) and (c) we show the amplitudes of the wave pack- 
ets ip±(q,t) corresponding to fig. 0] (a), while in fig. [5] 

(b) and (d) we display the same but for an initial Fock 
state with n = 2 photons. The Rabi model wave packets 
are displayed in (a) and (b), while (c) and (d) show the 
wave packets for the regular JC model. Note that in (d) 
the wave packet either has two or three peaks depending 
on if the field contains 1 or 2 photons respectively. The 
imperfect revivals come from the weak coupling between 
the two oscillators, which induces the level splitting of 
the oscillators at q = 0. 

Figure UJ gives the effect of the p dependent coupling of 
the two displaced oscillators in the JC model. It shows 
the amplitudes \tp±(q,t)\ (upper plots) and IV'T.J.G/i*)! 
(lower plots), for the Rabi model with initial Fock state 
|n = 4). The parameters are as in fig. Q] (d) and it is 
clear how the lower adiabatic wave packet slides down 
the displaced oscillators, while the upper adiabatic wave 
packet is squeezed at first. In the JC model, the wave 
packets would look very similar to those of fig. [5] (c) and 
(d), since we know that the population Rabi oscillates 
between the two Fock states \n — 4) and \n = 5). We 
can thus conclude that the spreading and squeezing of 
the wave packets in the Rabi model is due to the absence 
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Figure 5: The absolute amplitudes of the wave packets cor- 
responding to the previous fig. [4] (a) without RWA (a) and 
(b) and with RWA (c) and (d). In (a) and (c) the number of 
photons of the initial state is as in fig. [5] (a); n = 0, while in 
(b) and (d) n = 2. 



of the p dependent coupling. This has been verified nu- 
merically by artificially neglecting the p term from the JC 
model and a similar evolution as in the Rabi case is ob- 
tained. This effect of the p term obviously holds for any 
initial Fock state, also for large photon numbers n, which 
has a broad wave packet. This may seem counterintuitive 
since for large n the wave packet is non vanishing where 
the two energy curves are far apart and no coupling is 
expected to take part. However, as n increases, the wave 
packet width also increases in momentum space and the 
p-coupling between the curves becomes stronger. 



B. Collapse-revivals 

One of the more interesting phenomenon of the JC 
model is the collapse-revival effect Q. This is a direct 
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Figure 6: The four different absolute amplitudes for the Rabi 
model in the strong coupling regime corresponding to fig. Q] 
(a), go = 1 and = 0.1. 



consequence of the "graininess" of the field: Various parts 
of the state will Rabi oscillate with different frequencies, 
depending on the photon number n, leading to a collapse 
of various physical quantities, such as the atomic inver- 
sion (cr z ). However, after some particular time the parts 
may come back in phase and a revival occurs. Revivals 
for evolving wave packets exist due to the anharmonic- 
ity of the potential surface, and depending on the or- 
der n of the anharmonicity terms (~ x n ) different types 
of revivals are obtained, see the review (6{|. For initial 
squeezed sates, the JC model may also exhibit super and 
fractional revivals apart from the regular revivals [70| . 

In (33j it was shown, in the weak coupling limit, that 
the inclusion of the counter rotating terms to the dy- 
namics induces a fast oscillating pattern on top of the 
regular JC one. This is, however, in the weak coupling 
limit where we expect only small changes, while we ex- 
pect larger effects if the coupling is increased. In the 
papers (33l . [43| an approximate method was introduced 
to study the collapse-revivals in the Rabi model. In fig. 
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[7]the atomic inversion is plotted; (a) the Rabi model, (b) 
the adiabatic approximation and (c) the JC model. The 
parameters are in this case v = 4, ft = 5 and go = 0.3. 
Even though the fidelity (|36[) has a minimum as low as 
« 0.8, the agreement between the exact (a) and adiabatic 
approximated (b) results are very good. There is, how- 
ever, a great difference between the Rabi and JC results. 
The "plateaus" between revivals are not as pronounced 
in the Rabi model, and the revivals occur earlier. 

In the adiabatic regime the revivals may be easily un- 
derstood and one may estimate the revival time. The 
two adiabatic potentials (j3"2")l . for example in fig. Q] (a), 
are almost harmonic with frequencies slightly bigger and 
smaller than u>. As the initial coherent states evolve on 
these potentials they will first come out of phase and 
later return, giving the revival. The anharmonicity and 
the non adiabatic contributions lead to non perfect re- 
vivals. These arguments are confirmed in fig. [8j where 
the amplitudes of the wave packets for the Rabi model 
(a), and for the JC model in (b), both corresponding to 
fig. [71 are shown. A rough estimation of the revival time 
can be derived. In the adiabatic limit we may consider 
the wave packets as evolving on the disconnected adia- 
batic energy curves 



V. 



Rabi 



x 2 /ft 2 
= T ± VT 



2g 2 x 



2„2 



u'j^RabiX 2 ft , _^ 

2 2 [6 ° 



and the revival time given above is of dubious validity. 
We can conclude that the above analytical estimate is 
only valid in the very adiabatic limit, and as such, it is 
interesting; in the adiabatic limit the revival time does 
not depend on the initial state and consequently the av- 
erage number of photons n. 




and working with the boson operators for the JC model 



Figure 7: The atomic inversion for the Rabi model (a), adi- 
abatic approximation of the Rabi model (b) and JC model 
(c). Clearly the JC model does not describe the dynamics 
properly in this case. The parameters are O = 5, go — 0.3 
and v = 4. 



/Q2 q 
Hjc w Lua^a ± J — + g 2 aU sa uj±j C a*a ± — , (38) 

where we have expanded the square roots and introduced 
the renormalized frequencies 



U)± = UJ±,Bab% — <^±,JC 



lT f (39) 



The revival times then approximate 



tR 



2tt 



U>- 



(40) 



We see that this rough estimation gives identical revival 
times for both the JC and the Rabi models and it is in- 
dependent of the initial state. Using the parameters of 
fig. [7| we find tR w 174. In the Rabi case, instead of 
an analytically approximate result of the revival times 
given above, we numerically calculate the curvatures of 
the adiabatic potentials. Such an approach for the above 
example gives a revival time tn — 125.8, which is much 
better compared to fig. [7] However, it should be men- 
tioned that measuring the curvature is not absolute since 
it depends on the x-interval used. It is known that the re- 
vival time depends on the initial state of the field. This is 
easily understood, since the coupled harmonic oscillators 
depend on "momentum" p, which is directly related to 
the amplitude of the initial coherent state. For large pho- 
ton numbers, the adiabatic limit go\/n <C A breaks down, 



At some point, the adiabatic approximation breaks 
down and the above picture of the collapse-revivals fails. 
Still one may have interesting collapse-revival patters, as 
is seen in fig. [9l This figure shows the same as fig. [71 
but for the parameters ft = 0.2, v = 4 and go = 2. The 
revivals are present in the Rabi model and seem very sta- 
ble, while in the adiabatic and JC situations they die out 
fast. The amplitudes of the wave packets are displayed 
in fig. [TOl where (a) gives the Rabi case and (b) the JC 
situation. The plot clearly indicates the existence of the 
collapse-revival effect. As the adiabatic approximation is 
not valid we can no longer try to understand the evolu- 
tion from the an adiabatic point of view. However, since 
ft < go (n) one may see the wave packet as evolving along 
the diabatic energies described by the displaced harmonic 
oscillators. Now they share the same frequency u> (the 
anharmonicity comes from the avoided level crossing) , so 
that the oscillation periods are identical. Just like in fig. 
[5] (a), we see revivals in fig. [9] (a), which occur when 
the splitted wave packets return to their initial position. 
For other parameter regimes, where neither the adiabatic 
nor the diabatic frame holds, the inversion pattern for the 
Rabi model looks more random. 



C. Squeezing 

Another non-classical property of the field interact- 
ing with an atom is squeezing, see Q. The field is 
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Figure 8: The absolute amplitudes of the wave packets are 
shown in (a) corresponding to the Rabi case of the previous 
figure [3 and (b) for the JC model. The two parts first come 
out of phase (collapse) and then return back in phase (re- 
vival). It is seen that the revival occurs at later times with 
the JC, as seen also in the previous figure. 



said to be squeezed if one of the quadrature measures 
(Ag 2 ) = (q 2 ) - (q) 2 or (Ap 2 ) = (p 2 ) - (p) 2 is smaller 
than 1/2. The Squeezing in the Rabi model has been 
studied in 46], and it was found that squeezing for the 
vacuum is mostly pronounced when fl ^> go and g\ ^> CI. 
The second condition implies that the lower adiabatic en- 
ergy curve has two minima, as in fig. Q] (b)-(d). Dipolc 
squeezing for the Rabi model was considered in [51]. In 
the standard JC model, the amount of squeezing is im- 
proved for large amplitude initial states, meaning states 
with a large average number of photons. The asymp- 
totic solutions in these limits have been applied to study 
squeezing of the JC model [7l| . 

In fig. [11] the quadrature for q is shown for four dif- 
ferent cases. In the first two plots, (a) and (b), the state 
of the field is initially in vacuum and (a) go = 0.1 and 
fl = 0.2, while in (b) g = 0.1 and fl = 2. In |46|, it is 
shown that enhanced squeezing may be achieved if the 
atom is initially in a linear combination of its internal 
states, such that the wave packet is evolving mainly on 
one of the adiabatic energy curves. This comes from the 
fact that squeezing comes from the anharmonicity of the 
energy curves. Note that in both (a) and (b) the system 
is highly detuned (fl = 2, 0.2) resulting in a small squeez- 
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Figure 9: Same as fig. [7] but with fl — 0.2 and go = 2 instead. 
The revivals are most pronounced for the non-RWA case (a). 
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Figure 10: The absolute amplitudes of the original basis 
({|±}}) wave packets if)±(q,t) corresponding to the previous 
fig. [9] Plot (a) gives the amplitudes for the Rabi model and 
the (b) the same for the JC model. 



ing. In (c) and (d) the state of the field is initially in a 
coherent state with v = 4, and (c) go = 0.3 and = 5, 
and in (d) go = 0.2 and fi = 1. Thus, in (d) the system 
is in resonance, fl = uj. 

Squeezing also manifests itself in the various quasi 
phase distributions of the field. Here we give examples 
of the Q-function [72| . defined as Q(a) = 7r~ 1 (a|p/|a), 
where \a) is a coherent state, and pf = (+|p|+) + (— ) 
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Figure 11: The absolute quadrature {Ax 2 } for the JC (circles) 
and the Rabi (solid) models. In (a) go = 0.1, fl = 0.2 and 
v — 0, in (b) go = 0.1, Q = 2 and v = 0, in (c) go = 0.3, 
Q, = 5 and ^ = 4 and in (d) go = 0.2, = 1 and v = 4. 



is the field density operator and p the full atom-field den- 
sity operator. For a coherent state of the field, the Q- 
function is given by a symmetric Gaussian in the com- 
plex a-plane, while if it is squeezed, the Gaussian is no 
longer symmetric, but squeezed in some direction. The 
Q-function for various types of JC models has been stud- 
ied in [73[ and for the Rabi model in (34J. Figure [l"2l gives 
the Q-function at different times for the Rabi model (left 
plots) and for the JC model (right plots), with an initial 
coherent state with v = 4. In (a) and (b) go = 0.3 and 
Q, = 5, and therefor the system is largely detuned which 
is seen since the Q-function does not split into two clear 
parts, as in (c) and (d) where go = 1 and fi = 1. In the 
last two plots the Q-function clearly splits up, however 
it is seen that for the Rabi model the two parts of the 
function comes back in phase contrary to what is seen 
for the JC model. The total time evolved in (a) and (b) 
is tf = 400, while in (c) and (d) it is only tf — 6. Since 
the coupling is larger in (c) and (d) together with hav- 
ing the resonance condition, the characteristic timescales 
become much shorter. During the splitting the system 
quantities collapse and the revival occurs when the parts 
come back together again. Therefor, in the Rabi model a 
revival takes place very early, compared to the JC model. 
It should also be noted that at half the revival time in 
the Rabi model the field has split into a Schrodinger cat 
state [10j with approximately same phase but different 
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Figure 12: Evolution of the Q-function for the Rabi model 
(left plots) and the JC model (right plots). In (a) and (b) 
go = 0.3, Q — 5 and v = 4, while in (c) and (d) we display a 
resonant situation Q — 1, and with go — 1. Numbers indicate 
time ordering and the final times (indicated by up to 9 in (a) 
and (b) and by up to 5 in (c) and (d)) are tf = 400 for the 
upper plots and tf = 6 for the lower ones. 




Figure 13: The classically obtained phase space trajectories 
(|41[) corresponding to the previous example of fig. [12] (c) . 



amplitudes, which was discussed by Phoenix [3J]. This 
is not the case of the JC model, where the splitting is 
solely in phase and not amplitude. 

To understand this amplitude splitting we write a = 
^75 {ic ~ iPc) and approximate the real and imaginary 

parts of a by its classical variables; q c = and 



p c = — ^ ± . In this semi-classical approximation we use 
the adiabatic potential (f32|) as the classical potential gov- 

2 

erning the evolution H± = + V±(q c ), and taking into 
account the level splitting (replace q c by \q c \ in the last 
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Figure 14: The evolution of a Fock state \n = 6} in the Rabi 
model for different times and = 1 and go = 0.5. 



term), we find 



Pc 



(41) 



(Q2 + 2glqlY V^+Mtf 



where the =F-sign corresponds to which potential sheet 
V±(q c ) is being used. In fig. [13] we present the classi- 
cal phase space trajectories for the example of fig. [TH 
(c). The 'kinks' at q c = are due to the second term 
of V±(q) in eq. (|32|) . arising from the variations of the 
potential. For the JC system, on the other hand, the 
Hamilton equations of motion are of the form 



4c = f±(q c ,Pc), 
Pc = - f±{pc,q c ), 



(42) 



resulting in perfect circular trajectories due to the ro- 
tational symmetry in 2-d. Another way to explain this 
splitting is given in (47| , where one considers the energy 
manifolds 



P 



± 



9l 

4 



(43) 



The contours obtained by keeping e±{p,q) = £q, where 
£o is the initial energy, govern the adiabatic trajectories 
in phase space. 

As was noted above for the resonant situation and 
strong coupling, go ~ u>, the Rabi model gives very in- 
teresting behaviours, different from the JC model. This 
does not only happen for initial coherent field states. For 
a Fock state \n), the Q-function is a "ring" with radii 
\a\ = y/n and thus completely undetermined phase. In 
the JC model, a Fock state evolves according to Rabi os- 
cillations between \n) and \n + 1) (the atom is initially 



in the excited state |+)), and therefor, the Q-function 
will have a " breathing" amplitude motion, but the phase 
stays undetermined. For the Rabi model however, the 
phase may not stay equally distributed over 2-7T, and 
sort of cat states may form. This is verified in fig. Qj] 
where the Q-function is plotted for four different times, 
t = 0, 50, 62.5 and 75 for an initial Fock state with n = 6 
and Q = 1 and go = 0.5. Clearly, at some specific times t, 
the Q-function consists of four distinct blobs and at other 
times two. For other initial Fock states and other param- 
eters one may obtain Q-functions consisting of more than 
four blobs. These states could be of interest for quantum 
information processing [74j . 
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Figure 15: The entropy for an initial coherent state (a) as in 
fig. [7]but with v — 3, and a Fock state (b) with go = 1, fi = 1 
and n = 0. Circles display the entropy for the JC model and 
solid line the Rabi results. 



D. Entanglement 




Figure 16: The absolute value of the autocorrelation function 
for the JC (solid line) and the Rabi (dotted line) model, for 
an initial coherent state with v = 15 and go = 1 and Q = 1 
(a) and go = 0.3 and il — 5 (b). 



The entanglement shared between the field and the 
atom is well studied in the JC model [9|. For example, it 
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is known that for initial coherent states of the field, the 
atom is entangled with the field at all times except for 
half the revival times, where it disentangles from the field 



provided a large field amplitude [75J. Also for the Rabi 
model, entanglement has been analysed [3(|. There are 
many measures of entanglement between two interacting 
subsystems A and B, where one of the more established 
ones is the von Neumann entropy 



Sa = -Tr \ftA log (Pa)] 



(44) 



and pA = Trs [p] is the reduced density operator for sub- 
system A. For a pure disentangled state we have Sa = 0, 
and it can be shown that if the system starts out in a 
pure state Sa (t) = S*s(i) [Zf3]. For the two- level atom- 
field system studied here, defining the coefficients 



(45) 



c ™ (0 = / KC^ife*)^, 



we find the eigenvalues X± of the reduced density matrix 
for the atom as 



(46) 




a +Ei c «(*)i 2 (Ei c «(*)r- 1 ) 



and the entropy 



S A = -A + log(A + )-A_log(A_). 



(47) 



Other measures for entanglement between bipartite sys- 
tems are possible f?l\ . 

Examples of the entropy are presented in fig. [15] for an 
initial coherent state (a), corresponding to the example 
of fig [7] but with v — 3, and a Fock state (b) at resonance 
fi = 1 and go — 1. As expected, the entanglement is in 
general larger in the Rabi model. For the coherent states, 
we see the typical disentanglement at the revival time in 
both models. For a Fock state we know that in the JC 
model, the atom disentangle from the field periodically 
with the Rabi oscillations. This is not the case in the 
Rabi model, since a set of Fock states of the field will be 
coupled to the atom and the disentanglement must occur 
between all field states and the atom simultaneously. 



E. Autocorrelation function 



The autocorrelation function 
A(t, a )=<*(t)|*(*)>, 



t > s, 



(48) 



contains much information about the dynamics of a given 
system. The mathematical properties of the autocorre- 
lation function have been studied by numerous authors 




Figure 17: The absolute value squared of the fourier trans- 
formed autocorrelation functions of the previous fig. 1161 Dot- 
ted line corresponds to the Rabi model, while solid line the 
JC model. 



[781 ] . Here we will take time s = (and for simplicity de- 
note the autocorrelation function by A(t)), and thus, the 
autocorrelation function measures fidelity between the 
initial state and the time evolved state. Expanding the 
initial state in eigenstates of the Hamiltonian (assuming 
a discrete spectrum), \^ } = Sn "!^™), we get 



A d) = E 



(49) 



Fourier transforming the autocorrelation function then 
gives a direct relation to the spectrum, E n of the Hamil- 
tonian 



\A(e)\o,}2\ c n\ 2 S(e-E n ). 



(50) 



Note that in the numerical simulations, time t is given 
on a lattice and the lattice length L t and spacing dt de- 
termine the e-lattice according to the FFT, and also the 
widths of the (S-function peaks in the spectrum. 

In fig. [T5] we give an example of the autocorrelation 
function A(t) for an initial coherent state with ampli- 
tude v = 15, and parameters go = 1 and Q = 1 (a) and 
,9o = 0.3 and Q = 5 (b). The dotted line corresponds to 
the Rabi model while solid line to the JC model. We note 
that the peaks occur at around times t = 2mt (u) = 1), 
n = 0, 1, 2, and that there are two such reviving peaks 
for each 2mr and they become more separated for larger 
times. The autocorrelation function also contains the 
information about the revival times; when the distance 
between the two seperated peaks becomes equal to 2ir, 
interference of the wave packets take place resulting in a 
revival. Since the separation is faster for the JC model 
in the (a) plot, the revivals will occur earlier in the JC 
model compared to the Rabi model. The Rabi model 
seperation is not as smooth as for the JC model which 
will tend to smear out the revivals, which was already 
seen in fig. [7] In (b), corresponding to fig. \7\ the separa- 
tion speed is slightly larger for the Rabi model and the 
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revival times is therefor shorter than for the JC model. 
Measuring the revival times from fig. [TH] (a) and (b) for 
the JC model one finds tn « 94.2 and t R w 343.4 re- 
spectively, which coincide well with the known analytical 
approximate expression ._59| 



t R w 2ir- 



9o 



1 



1/2 



(51) 



Since the wave packet evolves on potential surfaces 
that are not fully harmonic, small spreading occur and 
the peaks in the autocorrelation function will as well 
broaden. This in turn leads to decreasing revival am- 
plitudes and the characteristic revival widths become 
longer, see fig. M ( c )- 

Figure [17] shows one part of the fourier transformed au- 
tocorrelation functions of fig. HH Clear peaks indicating 
the eigen enrgies are noted for both models. There seems 
to be two different energy spacings between the peaks, 
one slightly larger and one slightly lower than e = 1. In 
the JC model this is easily understood from the analytic 
expression (|3"0|) for the energies, which for large values 
on n can be linearly approximated over short n-intervals. 
From the fact that n 1 it follows that also q ^> 1 
and we are in the adiabatic regime, where the potential 
curves are well described by two centered harmonic os- 
cillators with frequencies slightly larger and smaller than 
one. This, of course, also applies to the Rabi model, and 
in [3^ | it is shown, in a perturbative way, that the spec- 
trum of the Rabi model approaches an equidistant one 
in the strong coupling limit. When n is decreased the 
linearization of the energies (|30|) is no longer justified, 
which as well follows from the fact that the adiabatic 
approximation breaks down in such regimes. 



VI. CONCLUSIONS 

In this paper we have used a wave packet approach to 
get a deeper understanding of two of the most commonly 
used systems in quantum optics and solid state physics; 
the JC and Rabi models. The wave packet method is 
purely numerical, but once the evolved wave packets are 
obtained, all various physical quantities are easily ob- 
tainable. Also, for the JC model, which is analytically 
solvable, one is often left with an infinite sum represent- 
ing some physical quantity, and either approximations 
or numerical methods are used to get a "closed" form 
of the quantity. More important, from the knowledge of 
how the wave packet evolves on the two coupled poten- 
tial curves, one can understand several phenomena such 
as Rabi oscillations, entanglement and collapse-revivals, 
and even find approximate results for some quantities in 



the adiabatic or diabatic limits. Another advantage of 
this approach is that it is numerically easily obtained for 
any parameter regimes, especially in the strong coupling 
regime, where other approaches such as truncation meth- 
ods tend to be cumbersome. 

The focus has been to study some of the known phe- 
nomena of the models and compare the two. Several 
differences have been observed, and, of course, especially 
when the RWA is known to break down. The comparison 
between the two models is, however, mainly to under- 
stand the wave packet evolutions, and not to point out 
the physical differences between the two. For very strong 
couplings, the Rabi model shows many interesting effects, 
for example splitting of the phase-space function into sev- 
eral parts, and also splitting in am plit ude. This effect has 
not been studied in detail earlier [34j. Another interest- 
ing aspect seen from the wave packet simulations is the 
effect of the "momentum" coupling between the two en- 
ergy curves, present in the JC model. It has been shown 
that this term bounds the wave packet close to the origin, 
making it not slide down the two energy curves. The link 
between the JC model and curve crossing problems is of 
great interest since knowledge from either of the two may 
give insight into the other one. The semi-classical treat- 
ment presented in subsection III CI shows how to derive 
time dependent two-level Hamiltonians, which have asp- 
plications in several fields of quantum mechanics [o2l 179 ] . 

The method of wave packet propagation may easily 
be applied to similar or extended models, for example 
Kerr systems [8 Of , pumped systems [81( , the Dicke model 
[721 . [HJ and multi- level models [83| and it is believed that 
much insight of these systems may be gained. The rela- 
tion between the models studied in this paper and models 
describing diatomic molecular dynamics is planed for a 
future project (see also 184]), and especially study the ef- 
fect of anharmonicity in the potential curves on the phys- 
ical quantities. For such analysis it is likely that shape 
invariant potential techniques (SUSY) may be applica- 
ble [HI]. Interference effects between the evolved wave 
packets, but for a molecular system, will be presented 
elsewhere [86| . These results may have interesting appli- 
cations in the JC or Rabi model. 
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